Digital FM demodulator using a lagrangian interpolation function

ABSTRACT

In a preferred embodiment, an FM demodulator includes an I and Q sampler for generating data samples of I and Q signals of an FM signal. A signal processor derives samples of a modulating waveform of the FM signal from the I and Q data samples, using a Lagrangian interpolation function of the I and Q data samples. Variables include derivatives of the I and Q signals obtained by differentiation of the Lagrangian function. Preferably, the processor further utilizes a correction factor to correct the modulating waveform samples obtained from the Lagrangian-based interpolation, to thereby derive a modulating waveform exhibiting low signal distortion.

FIELD OF THE INVENTION

The present invention relates generally to apparatus and methods for demodulation of a frequency-modulated (FM) signal, and in particular, to a computationally efficient, low distortion demodulation of a digitally sampled FM signal.

BACKGROUND OF THE INVENTION

Frequency modulation ("FM") is widely used in radio broadcasting as well as in a variety of other applications such as wireless telecommunications and musical recording. In FM systems, a sine-wave carrier is modulated so that its instantaneous frequency differs from the carrier frequency by an amount that is proportional to the instantaneous amplitude of the modulating wave. The modulated carrier is then demodulated at its destination.

A relatively modem method of demodulating an FM signal is generally referred to as digital FM demodulation. With this method, the FM signal to be demodulated is applied to an analog to digital (A/D) converter where it is sampled at a periodic rate. The sampled FM signal is then typically separated into time varying in-phase (I) and quadrature-phase (Q) signals, where the Q signal is the I signal delayed in phase by 90° at the carrier frequency. A digital signal processor is then employed to recover the modulating wave of the FM signal from the sampled I and Q signals.

Conventionally, the digital signal processor recovers the modulating wave by computing the function ##EQU1## From a computational standpoint, this approach has several disadvantages resulting in increased processing time and/or signal processor complexity. One drawback is that the evaluation of the arctangent function either requires a look-up table, or the use of an approximation algorithm. Moreover, differentiation of the arctangent function requires mapping of the output to the correct quadrant.

Accordingly, since it is desirable to keep signal processor complexity to a minimum in digital FM applications, there exists a need for computationally efficient processing of the time sampled I and Q signals to recover the information within the FM signal.

SUMMARY OF THE INVENTION

Embodiments of the present invention can provide computationally efficient, digital FM demodulation with low signal distortion. In a preferred embodiment, an FM demodulator includes an I and Q sampler for generating data samples of I and Q signals of an FM signal. A signal processor derives samples of a modulating waveform of the FM signal from the I and Q data samples, using a Lagrangian interpolation function of the I and Q data samples. Variables include derivatives of the I and Q signals obtained by differentiation of the Lagrangian function. Preferably, the processor further utilizes a correction factor to correct the modulating waveform samples obtained from the Lagrangian-based interpolation, to thereby derive a modulating waveform exhibiting low signal distortion.

BRIEF DESCRIPTION OF THE DRAWINGS

For a full understanding of the present invention, reference is had to exemplary embodiments thereof, considered in conjunction with the accompanying drawings like reference numerals depict like elements or features, for which:

FIG. 1 is a block diagram of an embodiment of an FM demodulator in accordance with the present invention;

FIGS. 2A and 2B show configurations for the I and Q sampler of FIG. 1;

FIGS. 3 and 4 show performance curves for a prior art FM demodulator;

FIGS. 5-11 show performance results for various embodiments of FM demodulators in accordance with the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 illustrates a block diagram of an exemplary embodiment of an FM demodulator in accordance with the present invention, designated generally as 10. Demodulator 10 may form a portion of a receiver in a wireless telecommunications system, e.g., as part of a base station receiver in a cellular system. It is understood that demodulator 10 may also be employed within other types of FM systems.

An input FM signal v(t) to be demodulated is applied to input line 11 of I & Q sampler 12. I & Q sampler 12 provides a sampled data in-phase (I) component of v(t) on line 13 and a sampled data quadrature-phase (Q) component of v(t) on line 14. The sampled data I and Q components are data streams . . . I_(n-1), I_(n), I_(n+1), . . . and . . . Q_(n-1), Q_(n), Q_(n+1) . . . , respectively, where successive data samples within each data stream such as I_(n), I_(n+1) are spaced from one another in time by T seconds.

Suitable circuit configurations for I & Q sampler 12 are well-known in the art. FIG. 2A shows an exemplary I & Q sampler, designated generally as 12a. The FM signal v(t) is sampled by A/D converter 22 to provide a sampled FM signal on line 23. Sampling is performed in a periodic manner and a sufficiently high sampling rate to avoid loss of signal information. The sampled FM signal is applied to each of first and second multipliers 24 and 26, respectively, which multiply the signal by respective functions cos(ω_(c) t) and -sin(ω_(c) t) in which ω_(c) is the carrier frequency (in rads⁻¹) of the FM signal v(t). The multiplied outputs are filtered by respective low pass filters 27 and 28 to provide the sampled data I and Q signals on output lines 13 and 14, respectively.

FIG. 2B depicts an alternative configuration for the I & Q sampler, designated generally as 12b. The analog FM signal v(t) is applied to each of multipliers 24 and 26, where the signal is multiplied by cos(ω_(c) t) and -sin(ω_(c) t) functions, respectively. The multiplied outputs are filtered by filters 27 and 28 and then sampled by respective A/D converters 29 and 31 to provide the sampled data I and Q signals.

Referring again to FIG. 1, the sampled data I and Q signals on respective lines 13 and 14 are applied to signal processor 15, which reconstructs the modulating wave . . . m((n-1)T), m(nT), m((n+1)T), . . . , where m(t) denotes the modulating wave of the FM signal v(t). The sampled m(nT) data stream is then applied to the digital-to-analog (D/A) converter 20, which provides the demodulated analog output signal m(t) on line 21. If desired, the output signal can be applied to a standard reconstruction filter, 18, to smooth out the signal. Optionally, memory, 19, can be incorporated in the signal processor to store the completely sampled modulating waveform for subsequent conversion by D/A convertor 20 when needed.

In wireless telecommunications and other applications, signal processor 15 also provides a sampled data signal on line 17, corresponding to (I² +Q²), which is indicative of the power level of the received FM signal and therefore provides a useful Received Signal Strength Indicator (RSSI).

To provide a foundation for the signal processing of the sampled data I & Q signals by signal processor 15 to obtain m(t), certain aspects of FM signals will first be presented. In general, an FM signal v(t) may be expressed as:

    v(t)=v.sub.m cos  ω.sub.c t+φ(t)!,               (1)

where v(t) is the RF voltage, V_(m) is the amplitude of the RF carrier, and ##EQU2## where m(t) is the modulating waveform or message signal and K_(P) is the peak frequency deviation in rads⁻¹ /volt. The RF bandwidth "BW" of this signal can be estimated using the well-known Carson's role, which states that:

    BW=2K.sub.p f.sub.b                                        (3)

where f_(b) is the highest frequency component present in the baseband message signal, m(t).

To demodulate a received FM signal v(t), one can first exponentially modulate v(t) to provide a complex signal u(t), defined by:

    u(t)=v(t)e.sup.-jω.sbsp.c.sup.t.                     (4)

The real and imaginary, or in-phase (I(t)) and quadri-phase (Q(t)) components, respectively, of the signal u(t), are:

    I(t)=v(t) cos  ω.sub.c t!;                           (5a)

    Q(t)=-v(t) sin  ω.sub.c t!;                          (5b)

Accordingly, by employing the multiplicands cos (ω_(c) t) and -sin (ω_(c) t) in FIGS. 2A or 2B, the I and Q signals are readily obtained. (It is noted that frequency error of the FM signal between the transmitter and receiver are being ignored in this discussion. A constant frequency error would simply map to a d.c. bias at the output, and would not distort the demodulated output).

In analog form, I(t) and Q(t) may be expressed as: ##EQU3##

With prior art FM demodulators, the message signal is recovered using: ##EQU4##

However, as explained previously, from a processing standpoint there are disadvantages associated with computing the derivative of an arctangent function.

An alternative expression for the modulating function m(t), also employed in prior art FM demodulators, is: ##EQU5##

See, for example, U.S. Pat. No. 4,910,469, issued to Takahashi, which discloses an FM demodulator generally using the approach of Eqn. (9).

From a signal processing standpoint, the Takahashi approach has the advantages of affording simple numerical operations and also, of avoiding function mapping.

When equation (9) is transformed into the digital domain, deviations from an ideal, zero distortion demodulation occur as a result of finite sample intervals. The implementation of the two differentiators dI/dt and dQ/dt is, therefore, a critical aspect of the digital demodulation process.

By using a first-order approximation to differentiation, one can obtain: ##EQU6## where T is the sampling interval between successive samples of I and Q, such as between samples I_(n) and I_(n-1). A digital solution "R₁ " of Eqn. (9) can then be obtained as: ##EQU7## A complete reconstruction of m(t) is realized by evaluating the variable R₁ for each successive pair of time samples of the I and Q signals.

Referring to FIG. 3, it is apparent that an FM signal demodulated in accordance with Eqn. (11) has an output frequency deviation K_(O) that accurately tracks the input frequency deviation K_(I) of the modulated FM signal only for relatively low values of K_(I). Curve 32 shows an ideal response in which K_(O) =K_(I) for all values of K_(I), while curve 34 is a plot of K_(O) vs. K_(I) resulting from the approximation of Eqn. (11). In this example, a 1 kHz audio tone is frequency modulated with a peak frequency deviation K_(P) of 12 kHz.

As illustrated in FIG. 4, signal compression at high frequency deviations causes signal distortion. Curve 36 illustrates an ideal demodulated audio tone. Curve 38 shows a distorted audio tone representing the demodulated signal, realized with the Eqn. (11) approach.

In accordance with a first embodiment of the present invention, the performance of the above differentiator is improved by further configuring signal processor 15 to perform the following series expansion for sin⁻¹ (χ): ##EQU8## where R₁ is known from Eqn. (11), n is an integer, and nT corresponds to the point in time in which the nth samples of I(t) and Q(t) are taken. The modulating function m(nT) then equals R_(1c) /K_(P).

The performance improvement using only the first, third and fifth order terms of Eqn. (12) is depicted in FIGS. 5 and 6. Curve 50 of FIG. 5 shows the improved performance using the approach of Eqn. (12); curve 60 of FIG. 6 illustrates the corresponding distortion improvement. The results for FIG. 6 are for a 1 kHz tone that is frequency modulated with a peak frequency deviation K_(P) of 12 kHz. It is understood that other values of K_(P) can be used, depending on the application. For example, K_(P) may be as high as 15 kHz for wireless applications with 30 kHz wide channels, or 10 kHz or lower for 25 kHz wide channels.

In accordance with an alternate embodiment of the present invention, the FM signal is demodulated with the use of Lagrangian interpolation and subsequent differentiation. With this approach, interpolation functions for Q(t) and I(t) are obtained and differentiated to yield respective derivatives Q'(t) and I'(t). These derivatives are then substituted into equation (9) to yield a more accurate representation of the modulating wave.

Lagrangian interpolation can be used to derive an expression for a function f(x) based upon knowledge of equally spaced samples f₀, f₁, . . . f_(p) of the function f(x).

The Lagrangian interpolation formula can be expressed as: ##EQU9## where

    l.sub.o (x)=(x-x.sub.1)(x-x.sub.2) . . . (x-x.sub.p)

    l.sub.k (x)=(x-x.sub.0) . . . (x-x.sub.k-1)(x-x.sub.k+1) . . . (x-x.sub.p)

    l.sub.n (x)=(x-x.sub.0 ((x-x.sub.1) . . . (x-x.sub.p-1).   (14)

Thus, for example, if a 2nd-order (m=2) Lagrangian interpolation is employed to approximate the I(t) and Q(t) components of the FM signal, then ##EQU10## where Q_(n), Q_(n-1) and Q_(n-2) are the nth, (n-1)st and (n-2)nd samples of Q(t) at times t_(n), t_(n-1) and t_(n-2), respectively, and t is a point in time between t_(n-2) and t_(n). An analogous expression can be obtained for I(t).

Differentiating Eqn. (15) and the analogous expression for I(t) to obtain derivatives dQ(t)/dt and dI(t)/dt, and substituting into Eqn. (9) yields the following expression for the modulating function m(nT): ##EQU11##

Thus, the signal processor can employ the algorithm of Eqn. (16) to derive time samples m(nT)=R₂ /K_(P) of the modulating wave m(t), by operating on three successive samples at a time of I(t)--I_(n), I_(n-1) and I_(n-2) --and three successive samples of Q(t) at a time--Q_(n), Q_(n-1) and Q_(n-2). A complete reconstruction of m(t) is realized by evaluating the variable R₂ for each successive tri-sample set--e.g., Q_(n), Q_(n-1), Q_(n-2) to obtain m(nT); Q_(n+1), Q_(n), Q_(n-1) to obtain m((n+1)T); and so on.

In any event, it is understood that higher order Lagrangian interpolation functions can alternatively be derived and differentiated to yield Q'(t) and I'(t), and that the invention is not limited to 2nd order Lagrangian interpolation. For example, if a 3rd-order interpolation were to be used, four consecutive samples of I(t) and four consecutive samples of Q(t) would be necessary to derive each sample of m(t).

FIG. 7 depicts the performance of a FM demodulator employing the algorithm of Eqn. (16). Curve 70 shows the actual response of the demodulator as compared to the ideal response shown by curve 32. The output frequency deviation K₀ accurately tracks the input frequency deviation K_(I) up to values of K_(I) of about 3 kHz. For higher values of K_(I), the output frequency deviation K_(O) exceeds the desired response, producing signal expansion. Although the relative amplitude error performance is similar to that of the uncorrected first order differentiator described above in reference to FIG. 3, the signal is no longer compressed with high frequency deviations, which results in less signal distortion.

FIG. 8 depicts the resulting signal distortion of a demodulator using the approach of Eqn. (16). Curve 80 shows the actual response of a demodulated FM signal originally modulated with a 1 kHz tone, using a peak frequency deviation K_(P) of 12 kHz, while curve 36 shows the ideal response. Signal distortion is significantly improved as compared to the uncorrected first order differentiator of FIGS. 3-4.

In accordance with another embodiment of the subject invention, the modulating wave m(t) is recovered by employing a correction factor to correct the solution R₂ of Eqn. (16). The correction factor is of the form:

    R.sub.2C =R.sub.2 (1-β|R.sub.2 |)   (17)

where β can be selected to minimize the error at a particular input frequency deviation K_(I), or over a selected range of input frequency deviations K_(I). Samples of the modulating wave m(t) are then obtained using

    m(nT)=R.sub.2C /K.sub.P.                                   (18)

FIGS. 9 and 10 illustrate the performance of an exemplary FM demodulator that utilizes Eqns. (17) and (18) in the signal processing used for the demodulation process. In FIG. 9, the actual response of output frequency deviation K₀ vs. input frequency deviation K_(I) represented by curve 90, is seen to closely track the ideal response of curve 32 up to values of K_(I) greater than 10 kHz. This represents substantial improvement in linearity for high peak frequency deviations as compared to prior art digital FM demodulators, such as that described in reference to FIG. 3. FIG. 10 shows the actual amplitude response, curve 100, using this demodulator for a 1 kHz tone (curve 36), where the 1 kHz tone was originally modulated using a peak frequency deviation K_(P) of 12 kHz. From the curves it is apparent that linearity is dramatically improved using this approach. It is noted that a β value of 2.46×10⁻⁶ was selected in the FM demodulator to obtain the performance depicted in FIGS. 9 and 10. This value was chosen to minimize error in the derived modulating waveform m(nT) at corresponding values of K_(I) from 5-10 kHz.

A comparison between the signal distortion of four FM demodulators, each employing a corresponding one of the four algorithms disclosed herein, is given in Table 1 below. The results are for a 1 kHz tone frequency modulated with a peak frequency deviation of 12 kHz, as discussed above. In terms of maximum amplitude error and mean square error, the corrected 2nd-order Lagrange differentiator approach (the R_(2C) solution) provides superior results. Accordingly, the most preferred embodiment employs the R_(2C) approach.

                  TABLE 1                                                          ______________________________________                                                    Maximum Amplitude                                                                            Mean Square Error                                     Algorithm  Error (dB)    with 1 kHz Tone                                       ______________________________________                                         R.sub.1    -2.67         0.0078                                                R.sub.1c   -0.94         0.0041                                                R.sub.2    +3.02         0.0017                                                R.sub.2c   -0.39         0.0009                                                ______________________________________                                    

FIG. 11 depicts the signal to distortion results for each of the four FM demodulators as a function of digital signal processor (DSP) instruction cycles, for a 1 kHz tone and K_(P) =12 kHz, as before. While the computation time for the demodulator employing the R_(2C) approach is approximately 41% longer than one employing the prior art R₁ approach, signal to distortion improves from about 18 dB to 27.5 dB. The demodulators using the R_(1C) and R₂ solutions also significantly improve distortion while trading off computation time.

It will be understood that the embodiments disclosed herein are merely exemplary and that one skilled in the art can make many modifications and variations to the disclosed embodiments without departing from the spirit and scope of the invention. All such modifications and variations are intended to be included within the scope of the invention as defined by the appended claims. 

What is claimed is:
 1. A method for demodulating a frequency-modulated (FM) signal, comprising:generating data samples of in-phase (I) and quadrature-phase (Q) components of the FM signal; and reconstructing the original modulating waveform of the FM signal from said I and Q components, using a Lagrangian interpolation function of said I and Q components.
 2. The method according to claim 1, wherein said modulating waveform is derived using the following equation: ##EQU12## where, n is an integer;I_(n) and Q_(n) are the nth data samples of the I and Q components of the FM signal, respectively; T is the sampling interval between successive I component samples and between successive Q component samples of the FM signal; I_(n) ^(') and Q_(n) ^(') are the derivatives of the I and Q components, respectively, of the FM signal, at sample n; K_(P) is a peak frequency deviation of said FM signal; and m(nT) is the derived sample of said modulating waveform at a time corresponding to nT seconds.
 3. The method according to claim 2, wherein said Lagrangian interpolation function is a second-order interpolation function, and wherein said time samples of said modulating waveform are derived using the following equation: ##EQU13## where, I_(n-1) and Q_(n-1) are the (n-1)st data samples of said I and Q components of said FM signal; andI_(n-2) and Q_(n-2) are the (n-2)nd data samples of said I and Q components of said FM signal.
 4. The method according to claim 3, further comprising correcting said samples m(nT) of said modulating waveform by means of the following equation:

    m(nT)=R.sub.2C /K.sub.P,

where

    R.sub.2C =R.sub.2 (1-β|R.sub.2 |),

where β is a constant selected to minimize error in said modulating waveform.
 5. The method according to claim 2, wherein said peak frequency deviation is in the range of about 10-15 kHz.
 6. The method according to claim 1, wherein said generating data samples of I and Q components of said FM signal includes:sampling said FM signal to provide a sampled data FM signal; and generating said data samples of said I and Q components from said sampled data FM signal.
 7. The method according to claim 1, wherein said generating data samples of I and Q components of said FM signal comprises:generating analog I and Q components of said FM signal; and sampling said analog I and Q components to generate said data samples of said I and Q components of said FM signal.
 8. The method according to claim 2, wherein said FM signal is a received FM signal in a wireless telecommunications system.
 9. The method according to claim 8, further comprising the step of outputting of a term (I_(n) ² +Q_(n) ²) as a received signal strength indicator (RSSI).
 10. A method for demodulating a frequency modulated (FM) signal, comprising:generating data samples of in-phase (I) and quadrature-phase (Q) components of the FM signal; recontructing the original modulating waveform of the FM signal from said I and Q components, using the following equation:

    K.sub.P m(nT)=R.sub.1 +R.sub.1.sup.3 /6+3R.sub.1.sup.5 /40+. . . ,

where, n is an integer; T is the sampling interval between successive I component samples and between successive Q component samples of the FM signal; m(nT) is the derived sample of said modulating waveform at a time corresponding to nT seconds; R₁ is defined by the equation: ##EQU14## where I_(n) and Q_(n) are the nth time samples of the I and Q components of the FM signal, respectively; and I_(n-1) and Q_(n-1) are the (n-1)st time samples of the I and Q components of the FM signal, respectively.
 11. An apparatus for demodulating a frequency modulated (FM) signal, comprising:an I and Q sampler for generating data samples of an in-phase (I) and a quadrature-phase (Q) component of the FM signal; and a signal processor for reconstructing the original modulating waveform of the FM signal from said I and Q components, using a Lagrangian interpolation function of said I and Q components.
 12. The apparatus according to claim 11, wherein said signal processor is operative to reconstruct the modulating waveform using the following equation: ##EQU15## where, n is an integer;I_(n) and Q_(n) are the nth data samples of the I and Q components of the FM signal, respectively; T is the sampling interval between successive I component samples and between successive Q component samples of the FM signal; I_(n) ^(') and Q_(n) ^(') are the derivatives of the I and Q components, respectively of the FM signal, at a time corresponding to nT seconds, and based on differentiation of said Lagrangian interpolation function of respective of said I and Q data samples; K_(P) is a peak frequency deviation of said FM signal; and m(nT) is the derived sample of said modulating waveform at a time corresponding to nT seconds.
 13. The apparatus according to claim 12, wherein said Lagrangian interpolation function is a second-order interpolation function, and wherein said signal processor is operable to derive said samples of said modulating waveform using the following equation: ##EQU16## where, I_(n-1) and Q_(n-1) are the (n-1)st data samples of said I and Q components of said FM signal; andI_(n-2) and Q_(n-2) are the (n-2)nd data samples of said I and Q components of said FM signal.
 14. The apparatus according to claim 13, wherein said signal processor is further operable to correct said samples m(nT) of said modulating waveform by means of the following equation:

    m(nT)=R.sub.2C /K.sub.P,

where

    R.sub.2C =R.sub.2 (1-β|R.sub.2 |),

where β is a constant selected to minimize error in said modulating waveform.
 15. The apparatus according to claim 11, wherein said I and Q sampler comprises:an analog to digital (A/D) converter for sampling the FM signal to provide a sampled FM signal; first and second multipliers arranged in a parallel circuit configuration for multiplying the sampled FM signal by respective of cos(ω_(c) t) and -sin(ω_(c) t) functions to provide first and second multiplied signals, respectively, where ω_(c) is a carrier frequency of the FM signal and t is time; and first and second low pass filters for filtering the respective first and second multiplied signals to provide said I and Q data samples, respectively.
 16. The apparatus according to claim 11, wherein said I and Q sampler comprises:first and second multipliers arranged in a parallel circuit configuration for multiplying the FM signal by respective of a cos(ω_(c) t) and -sin(ω_(c) t) function to provide first and second multiplied signals, respectively, where ω_(c) is the carrier frequency of the FM signal and t is time; first and second low pass filters for filtering the respective first and second multiplied signals to provide first and second filtered signals, respectively; and first and second analog to digital (A/D) converters for sampling the respective first and second filtered signals to provide respective of said I and Q data samples.
 17. The apparatus according to claim 11, wherein said signal processor is operable to provide an output signal indicative of magnitude of said FM signal.
 18. The apparatus according to claim 11, further including a digital to analog (D/A) converter coupled to said signal processor for converting said samples of said modulating waveform into an analog modulating waveform of said FM signal.
 19. The apparatus according to claim 18, wherein the signal processor includes a memory, said memory operable to store a complete sampled modulating waveform for subsequent D/A conversion by said D/A converter.
 20. An apparatus for demodulating a frequency modulated (FM) signal, comprising:circuit means for generating samples of in-phase (I) and quadrature phase (Q) components of the FM signal; and processing means for reconstructing the original modulating waveform of the FM signal from said I and Q components, using a Lagrangian interpolation function of said I and Q components. 